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Express Mail Label No. EV 320890794 US 

[0001] EXTENDED ALGORITHM DATA ESTIMATOR 

[0002] CROSS REFERENCE TO RELATED APPLICATION(S) 

[0003] This application claims priority from U.S. Provisional Application No. 
60/409,973, filed on September 9, 2002, which is incorporated by reference as if fully 
set forth. 

[0004] FIELD OF INVENTION 

[0005] The invention generally relates to data estimation for wireless 
communication systems. More particularly, the invention relates to an extended 
algorithm (EA) for data estimation in such systems. 

[0006] BACKGROUND 

[0007] In some proposed wireless systems, data is wirelessly transmitted block 
by block with a separation interval between successive blocks. This property permits 
the application of joint detection (JD) in receivers to suppress the inter-symbol 
interference (ISI) and multiple-access interference (MAI). Single user detectors (SUDs) 
are used to estimate data of signals that go through a single downlink channel. An 
advantage of the SUD is that it can be implemented eflSciently by fast Fourier 
transform (FFT), based on the rationale that a square Toeplitz matrix can be 
approximated as a circulant counterpart of the same size. 

[0008] When Toeplitz matrices are shortened along their longer dimensions to 

square matrices and replaced with their circulant counterparts, an approximation error 
is introduced. This error is most prevalent in the head and tail portions of the matrix. 
In many systems, the data associated with the head and tail portions carries system 
information required by the receivers, such as power control bits and the transport 
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format combination indicator (TFCI) in the proposed third generation partnership 
project (3GPP) wideband code division multiple access (WCDMA) time division duplex 
(TDD) system. It is desirable to enhance data estimation in such systems. 

[0009] SUMMARY 

[0010] The present invention provides a computationally efficient and accurate 
implementation of a data estimator for systems, such as Frequency Division Duplex 
(FDD) or Time Division Duplex (TDD) Code Division Multiple Access (CDIVIA) systems. 
Described herein is a method of implementing a data estimator wherein no data of 
interest for estimation will be affected significantly by the circular approximation 
error. To achieve this, all of the square circulant matrices are extended. The 
advantages of the extended approach arise fi-om two aspects: (1) the avoidance of loss of 
multipath signals in the tail part of each data field and (2) the avoidance of error due to 
the Toeplitz to circulant matrix transformation. As a result, longer discrete Fourier 
Transforms (DFTs) or FFTs are performed when implementing the EA. In order to 
minimize the required computations with DFT, the extended sizes should preferably be 
limited djoiamically to their lower bounds according to specific data block length and 
channel delay spread. However, if prime factor algorithm (PFA) is used, increasing the 
FFT length usually does not result in the increase of computational complexity. It is to 
be noted that the computations can be minimized through choosing a proper FFT 
length in a certain range. In this case, fixed single-length FFTs are desirable by 
considering the longest block length and delay spread. The single-length FFT with 
PFA makes it possible that different lengths of data block (burst types) are supported 
by one algorithm only. This simplifies the implementation further since one piece of 
hardware is needed to cope with the single algorithm. 
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[0011] BRIEF DESCRIPTION OF THE DRAWING(S) 

[0012] A more detailed understanding of tlie invention may be had from tlie 

following description of preferred embodiments given by way of example only and to be 

understood with reference to the accompanpng drawing wherein: 

[0013] Figure 1 is a block diagram of a system used for implementing an EA with 

over-sampling in accordance with a preferred embodiment of the present invention; 

[0014] Figures 2A, 2B and 2C, taken together, are a flow chart illustrating the 

method steps implemented by the EA of Figure 1. 

[0015] Figure 3 is an illustration of raw block error rate (BER) of TFCI-1 versus 
signal to noise (SNR) in Case 1 channel (TFCI represents transport format combination 
indicator); 

[0016] Figure 4 is an illustration of raw BER of TFCI-2 versus signal to noise 
ratio (SNR) per code in Case 1 channel; 

[0017] Figure 5 is an illustration of raw. BER of all bits versus SNR per code in 
Case 1 channel; 

[0018] Figure 6 is an illustration of raw BER of TFCI-1 versus SNR per code in 
Case-2 channel; 

[0019] Figure 7 is an illustration of raw BER of TFCI-2 versus SNR per code in 
Case-2 channel; 

[0020] Figure 8 is an illustration of all bits versus SNR per code in Case-2 
channel; and 

[0021] Figures 9A and 9B are receiver implementations using extended 

algorithm data detection, 

[0022] DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 

[0023] The invention is generally applicable for data estimation in CDMA 

systems, such as 3GPP TDD mode and time division synchronous code division 
multiple access (TD-SCDMA). However, the following description, as an example. 
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refers to a TDD CDMA system model and algorithm. In the described example, if a 
composite spread signal is sent from one transmitter site to one reception site, the 
received signal, r , is composite spread signal, 5 , passed through a single channel, H . 

is the channel response matrix. This process can be represented asr = ^ + n, 
where n is the noise vector. With W being the length of channel response, H_ takes the 
form per Equation 1. 

K 



Equation 1 



H is of size (L+W-1) by L. L denotes the number of chips in the time period of interest, 
such as a data field (block). The composite spread signal s can be expressed as ^ = Crf , 
where the symbol vector and the code matrix C take the form per Equation 2. 

d = (rfj ^"-y^KN y Equation 2 



T denotes the transposition and C is per Equation 3. 
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c = ]g!\c''\...,c'''\ 



Equation 3 



Each C^*^ is per Equation 4. 



Equation 4 



[0024] Q, K, and (=L/Q) denote the spread factor (SF), the number of active 

codes, and the number of sjnnbols carried on each channehzation code, respectively. 
[0025] The SUD comprises two stages: (a) channel equalization and (b) de- 
spreading. In the first stage, the composite spread signal s is estimated from 
r = HjS + n , preferably through a minimum mean squared error (MMSE) equalizer or a 
zero forcing solution. 
A MMSE equalizer is per Equation 5. 



1 = [Rh -^^^iYh% 



Equation 5 



A zero forcing solution is per Equation 6. 
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I = Rh'^U." r , Equation 6 

I is the identity matrix, i?„ = H_" H_issi square Toeplitz matrix of size L, per Equation 
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Equation 7 

' * ' denoting the conjugate operation. In the second stage, a simple de-spreading 

process is performed to estimate the symbol sequence , , per Equation 8. 

d = C^s. Equations 

[0026] To implement Equation 5 efficiently, it is desirable to approximate the 
algorithm properly. To do this, first the Toeplitz matrix is extended from the size of 
L+W-1 by L to the size of Lm+W-1 by Lm and then the square matrix R„ is extended 
from the size of L to Lm with Lm > L+W-1, while keeping the banded and Toeplitz 
structure of the matrices intact. The vector r is extended to the length Lm by zero 
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padding, if the length of r is less than Lm. The vectors s and n are effectively 
automatically extended, due to the padding extension of the other vector/matrix. The 
extended versions of the matrices and vector , 5, r and n are denoted as H^^y 

^£ » §.E9 Le He 9 respectively. R^is as follows, = H"eKe • 
[00271 The last Lm-L elements of can all be deemed to be zero, which is 
essential in understanding the avoidance of the implementation error. With these 
notations, Equation 5 can be rewritten as Equation 9. 



is as follows =KeIe +5e- Because Equation 9 is the extended version of 

Equation 5, there should not be any differences between them if only the first L 
elements of are considered. The last W-1 rows of the H^^ are cut to get a new 
square matrix of size Lm, denoted by H ^ . Suppose R^.^ and H^cir represent the circular 
counterparts of R^ and H^^ , respectively. From R^ and , R^^^ and H^^ir 
constructed as per Equations 10 and 11. 




Equation 9 
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Equation 10 



iLcir 



K-i 



ho 



Equation 11 



Rem and ''^th of them are square matrices of size Lm. Because the last Lm-L 
elements of are zero, = Ks 1e He • Replacing , H_" and by /?c,>> Kar 



KsIe + 5e > respectively, in Equation 9, results in Equation 12. 
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r r 1-1 Equation 12 

Assuming H_s -Kor "^K^ (IL^ is an error matrix), Equation 12 can be expressed as 
Equation 13. 

(- . Equation 13 

^ = H_ciriLA^E is a column vector of length Lm and the error matrix is per Equation 
14. 



K- 



Equation 14 



The non-zero elements are located in a triangular area between the first W-1 rows and 
the last W-1 columns. Without the third term, Equation 13 is very similar to Equation 
5 in function when the first L elements are considered. The vector y is evaluated. 

According to the structures of the matrices H^iv ILt. » matrix X = HarKu, is a 
square matrix of size Lm with the structure per Equation 15. 
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X 



X X 



X 



X 



Equation 15 



X 



X X 



X 



X 



X 



X X 



The non-zero elements, denoted by 'x', are located only in two areas: (1) a triangular 
one between the first W-1 rows and the last W-1 columns and (2) a square one between 
the last W-1 rows and columns. Because the last (Lm-L) elements of are zero, all 

elements of y = H_cirKAlE ~ K§.e are zeroif Lm>L+W-l. When Lm>L+ W-1, Equation 
16 results. 



Equation 16 is a good approximation of Equation 5, when taking only the first L 
estimates into account. The first portion of Equation 16 is referred to as the EA. 
Similarly, the head and tail parts of a data field are affected significantly by replacing 
Toeplitz matrices with their circulant counterparts in Equation 5 directly without 
matrix extension. The implementation algorithm without matrix extension is called 
the truncated algorithm (TA). With TA, Equations 9 to 15 are still valid, except that 
Lm=L. The reason consists of two aspects: First, when replacing in Equation 5 with 




Equation 16 
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KcirOf size L, the length of the received signal vector r is preferably limited to be L. 
This results in the loss of multipath signals of the data in the tail part of a data field. 
Therefore, the estimates of the affected data become very poor. Second, when Lm=L, y 

is a column vector of length L, in which the first and last W-1 elements are non-zero 
elements. Because B^.^ =Rcir -^^^l is of banded structure along the diagonal, the 
inverse of B^.^ has approximately the same structure. Therefore, the relative large 
values of the column vector z = y are located in the first and last W-1 colunms and 
affect the estimates of | both in the head and tail areas. For the second reason, the 
number of the affected estimates depends on the channel response length W. The 
larger the delay spread of a channel (W) is, the more the estimates are affected. 
[0028] Additionally, the implementation algorithm given by Equation 16 can be 
extended to support over-sampling. With over-sampling, the effect of timing error is 
mitigated. Assuming sampling rate is M times of chip rate, M received signal vectors, 
denoted by r for m=l,2, . . . ,M, are available. However, the time interval between two 
successive samples in each r;^^ is the chip duration. Similarly, there are also M sets of 
channel response, denoted by h^"^^ = (\m A,m'""^vK-i,m) m=l,2,...,M. With these 
channel responses, a total of 2Mcirculant matrices Kor^m ^c/r.m for m=l,2,..., Mean 

be constructed. Accordingly, the implementation algorithm with over-sampling can be 
written as per Equation 17. 



m=:l 



T,H"ar^rf' Equation 17 



gI is the noise variance corresponding to the mth input vector . 

[0029] Before implementation, the value of Lm is determined. Because Lm is 
larger than L+W-1, Lm is chosen as per Equation 18. 

Lm = max {L}+max {W}+e, Equation 18 



-11- 



1-2-0394, lUS 



max {.} is the maximum value of {.} and 8 is a positive integer used to let Lm a good 
length for FFT implementation. For example, in UTRA wideband TDD system 
(WTDD), max {L}=1104, max {W}=114. e is chosen to be equal to 14 to let Lm=1232. 
With this length, FFT can be performed very efficiently by PFA because 1232 can be 
factored as 1232=7 x 1 1 x 16. With complex input, the real multiplies and adds required 
by the 1232-point FFT are 8836 and 44228, respectively. From Equation 18, Lm 
depends on specific system design. However, the implementation approach described 
herein is applicable to any other TDD systems such as UTRA narrowband TDD system 
(TD-SCDMA). 

[0030] In the following, the preferred implementation procedure of Equation 17 

is described as method steps under the assumption that the FFT length P equals the 
selected Lm. 

[0031] The first column g of the circulant matrix V R^,^ ^ + — cr^/ is computed 

based on the estimated channel response and noise power, 3delding Equation 19. 

Equation 19 



( 1 



Af / 1 

[0032] The circulant matrix ^ j R^.^ ^ ^ a7 domain is decomposed , 



yielding Equation 20. 



U^cir^m +^or^^ ] = ^XD^, Equation 20 



Dp and Dp are the P-point FFT and inverse FFT (IFFT) matrices defined as per 
Equation 21. 



DpX = ^xin)e~' >" and D;'x = -^x(n)e' >" , for k=0,l,...,P-l, 

Equation 21 
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is a diagonal matrix of size P, whose diagonal is Dpg. A^is denoted as 
A/e =diagiDpg). The relation between Dp and Dp is Dp^ ^{l/P)D*p. 
[0033] The circulant matrix Har^^ is decomposed in the FFT domain, 5delding 

Equation 22. 

Kcinm = D~p^H„ Dp y Equation 22 

A^^ is a diagonal matrix of size P, whose diagonal is DpU^ with 
Mm = Wm » Km » ' * ' ' K-\.m * ' ^of being the first column of Har,^ . 

[0034] The received signal vector r ^'"^ is reconstructed by zero padding to get the 

extended signal vector r^^^ of length P, 

[0035] The composite spread signal vector ? ^ is computed yielding Equation 23 
or Equation 24 in the frequency domain. 



f M 



' M M 

2iLar^r.f = D;'A-l2^\D,/r' Equation 23 



m=l m=l 



DpS_E = 



Y^iDpuJ ®(Dpr[^^^^^ and 1^ = D'p'iDpl,} , Equation24 

The operators ® and / denote the vector multiplication and division performed on 
element-by-element basis, respectively. The last P-L elements of ? ^ are rounded off to 
get another vector | of length L. 

[0036] The composite spread signal | is despread, 3delding?^ . 
[0037] Figure 1 is a block diagram of a system 100. For an over-sampled 

system, M sampled sequences are processed, . . . r^^^ and h^^^ .,.h^^\ For a chip rate 
sampled sequence, only one sampled sequence is processed, r^^^ and h^^^ . The system 
100 receives the signals . . . r at inputs 105i. . . 105m (105) and receives the channel 
impulse response h}^^ ... ^^^^ at inputs llOi . . .110m (110). The received signals r^*^ . . . 
r^^^ are zero padded in the tail by zero padding devices 115i...115m (115) until the 
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length of each sequence achieves length Lm. The extended sequences after zero 
padding are denoted as rf ... r^/^^ which exits block 115 via outputs 120i...120m (120). 
The channel impulse responses h}^\..h^^^ are zero padded in the tail by zero padding 
devices 125i...125m (125) until the length of the extended sequence achieves length 
Lm. The extended sequences after zero padding are denoted as Mi . . . « ^ which exits the 

zero padding devices 125 via outputs 130i... 130m (130). DFT or FFT blocks 135i...135m 
(135) receive the outputs 120 from zero padding devices 115 and perform DFT or FFT 
on /J^..r^/^ Fi/^)...F(/^^). DFT or FFT blocks 140i...140m (140) receive the 
outputs 130 from zero padding devices 125 and perform DFT or FFT on ...w^, 
F(mi)...F(m^). Conjugate devices 145i...145m (145) conjugate F(wi) ... F(wj^ ), 
F(mi)* ...F(m^)*. Element-to-element multipliers 150i...150m (150) multiply the 
sequences Firf),..F(/^') and F(u,)\-- F(u^)\ F{rf)^ F{uS F(rf')' F(u^,r , 
[0038] All of the M sampled sequence results are added element-to-element by 

M 

adder 175 with 5^F(r^"'^) F(w^)* , M = 1,2,. ..,M. A channel correlation vector £ is 
generated by a channel correlation generator 180 using extended channel response 

M 

sequences , with g =^g^'"^ . 

[0039] Using a IVIIMSE algorithm, a noise variance al is added to the first 

element of vector g^'"^ Vector g^'"^ is generated using . The i-th element of the 
vectors ^^"^ for the m-th sampled sequence is computed by: circulating the conjugate 

vector by downshifting i-1 elements and multiplsdng the shifted vector by the 
vector such that £^'"^(/) = w^.o^w^^m . m = l,2,...,Af . A DFT or FFT device 185 
performs a DFT or FFT on channel correlation vector g , F(g) . Divider 190 divides 
element-by-element the output of adder 175 by the output of DFT or FFT device 185, 
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such that — . Inverse DFT or inverse FFT device 194 is performed on 

XF(rr)-F(«„r 

the output of divider 190, such that F \— ) . The output of the inverse 

Fig) 

DFT or inverse FFT device 194 is an estimate of the composite spread signal s, 
Despreader 198 despreads the output of inverse DFT or inverse FFT device 194 to 
obtain the estimated data S3anbols d . 

[0040] Referring to Figures 2A, 2B and 2C, the procedure for performing an EA 
in accordance with a preferred embodiment of the present invention is described as 

follows: 

[0041] In step 205, system 100 receives signal r^^^ at input 105 and receives the 

channel impulse response ^^^^at input 110. 

[0042] In step 210, the received signal r^^^ is zero padded in the tail by zero 

padding device 115 until the length of sequence achieves length Lm. The extended 
sequence after zero padding is denoted as r^^^ which exits block 115 via output 120. 

[0043] In step 215, the channel impulse response h}^^ is zero padded in the tail 

by zero padding device 125 until the length of the extended sequence achieves length 
Lm. The extended sequence after zero padding is denoted as which exits the zero 
padding device 125 via output 130. 

[0044] In step 220, DFT or FFT block 135 receives the output 120 from zero 

padding device 115 and performs DFT or FFT on rf such that Firf) . DFT or FFT 
block 140 receives the output 130 from zero padding device 125 and performs DFT or 
FFT on w, , F(wj). 

[0045] In step 225, conjugate device 145 conjugates F(u^) , F(u^)* . 
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[0046] In step 230, element-to-element multiplier 150 multiplies the sequences 

F(rf) and F(w,)* producing F(rf)'F{u,)\ 

[0047] In step 235, for over-sampling system with M sampled sequences, steps 

210 to 230 are repeated for sampled sequences 2, . . . ,]VI, F(rf^ ) * F(u„ )* , for m = 2, . . . , IM. 
[0048] In step 240, all of the M sampled sequence results obtained in steps 230 

M 

and 235 are added element-to-element by adder 175, ^Fi/^^)- F(u^y , for m = 
1,2,.. .,M. 

[0049] In step 245, a channel correlation vector g is generated by a channel 

correlation generator 180 using extended channel response sequences u^,^"yUM y such 

[0050] In step 250, a DFT or FFT 185 performs DFT or FFT on channel 

correlation vector g , F(g) . 

[0051] In step 255, divider 190 divides element-by-element the result in step 240 

M 

J,Firf)-F(uJ 



by the result in step 250, 



m=l 



[0052] In step 260, an inverse DFT or inverse FFT 194 is performed on the 



result of step 255, F ^ 



producing the estimated composite spread 



signal, 5. 

[0053] In step 265, despreader 198 despreads the output of step 260 to obtain 

A 

the estimated data sjnnbols d . 
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[0054] In simulation, the described model is based on K = 12, and that codes are 

transmitted with equal code power and the effects of midamble signal on the second 
data field have been cancelled completely. Each code has a SF =16. A total 1104 chips 
in a data field are assumed (burst type 2 in WTDD). Because there are two data fields 
in a timeslot in WTDD, the last 8 bits (4 complex sjnnbols) of the first data field, and 
the first 8 bits of the second data field are defined as TFCI-1 and TFCI-2, respectively. 
Two algorithms, tnmcated algorithm (TA) and extended algorithm (EA), are employed.. 
The raw bit error rates of TFCI-1, TFCI-2 as well as all bits are evaluated for both EA 
and TA with chip rate sampling. 1000 timeslots are accumulated for each SNR point. 
The simulations are run over the WG4 Case-1 and 2 channels. 

[0055] Figures 3 and 4 present the performance of TFCI-1 and TFCI-2 in WG4 

Case-1 channel when EA and TA are used. A significant performance gap between EA 
and TA is found as shown in Figure 3, Because TFCI-1 locates in the tail part of the 
first data field, the performance degradation of the TFCI-1 with TA results from two 
reasons: (1) loss of TFCI-1 multipath signals and (2) Toeplitz to circulant matrix 
replacement in the implementation because the channel response length W is small 
(W=4) for WG4 Case-1 channel. This conclusion is confirmed by the results shown in 
Figure 4. Because the performance in Figure 4 is for the TFCI-2 which locates in the 
head part of the second data field, the most probable reason that affects its 
performance with TA must be the second one: matrix replacement. From Figure 4, it 
can be seen that the performance of TFCI-2 with EA and TA is nearly identical. This 
implies that the estimate error introduced to the TA through the matrix replacement is 
very limited because of the small value of W for WG4 Case-1 channel. For example, 
when W=4, only a quarter of the first symbol is affected because of the relationship 
SF=16. 

[0056] Figure 5 shows the raw BER of all bits when EA and TA are assumed. 

Comparing with EA, the major contributor to the loss of raw BER of all bits with TA is 
the TFCI-1 in each timeslot. 
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[0057] Figures 6 and 7 show the performance of TFCM and TFCI-2 in WG4 
Case-2 channel when EA and TA are adopted. Case-2 channel differs from Case-1 
channel by much larger delay spread (W=46) as well as stronger power of multipath 
signals. From Figure 6, the TFCI-1 with TA is seen to be almost destroyed both due to 
the loss of its multipath signals and matrix replacement. In Figure 7, the performance 
of TFCI-2 with TA is still much poorer than that with EA, which is due to the matrix 
replacement only. In this case, W=46 and hence the first threie s5nnbols (six bits) will 
be affected significantly. Figure 8 presents the raw BER of all bits for EA and TA in 
Case-2 channel. 

[0058] When a truncated algorithm (TA) is used in the implementation, the data 

in the head and tail parts of a data field are affected significantly by two aspects: the 
lost information of multipath signals due to cutting the channel response matrix to be 
square and the error due to replacing Toeplitz matrices with circulant ones. To 
overcome this problem, an EA is used. The EA avoids implementation errors by 
choosing the size of the extended matrices properly. To implement the EA with DFT, a 
djmamic-length EA is desirable, while for an EA with FFT of PFA, a fixed-length 
approach is more appropriate. In the fixed-length approach, the computation 
complexity can be minimized through choosing a proper PFA length in a certain range. 
The fixed-length EA makes different data blocls: lengths (burst types) be supported by 
one algorithm only. The fixed length EA simplifies the implementation further since 
one piece of hardware is needed to cope with the single algorithm. Simulation results 
show that the performance of the EA is much better than that of the TA, especially for 
the data in the head and tail parts of data fields. 

[0059] The invention can be implemented at a base station or wireless 
transmit/receive unit (WTRU), Hereafter, a wireless transmit/receive unit (WTRU) 
includes but is not limited to a user equipment, mobile station, fixed or mobile 
subscriber unit, pager, or any other t3^e of device capable of operating in a wireless 
environment. When referred to hereafter, a base station includes but is not limited to 
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a base station, Node-B, site controller, access point or other interfacing device in a 
wireless environment. 

[0060] Figures 9A and 9B are receiver implementations using extended 
algorithm data detection. Referring to Figure 9 A, radio frequency (RF) signals are 
received by an antenna 300. A sampling device 305 produces a chip rate received 
vector r. A channel estimation device 325 determines a channel impulse response h for 
the received vector. A single user detection device 310 uses the received vector r and 
the channel impulse response h to estimate the data vector d using the extended 
algorithm. The received vector r is processed by a channel equalizer 315 using the 
channel impulse response h to determine a spread vector s. A despreader 320 using the 
transmission codes C despreads the spread vector s to estimate the data vector d. 
[0061] Referring to Figure 9B, RF signals are received by an antenna 300. A 
sampling device 305 samples the received signal at a multiple M of the chip rate, 
producing M received vector sequences D...rm. A channel estimation device 325 
determines a channel impulse response 4i--.4m corresponding to each received vector 
Q. . .rm. A single user detection device 310 uses the received vector sequences q. . .rmand 
the channel impulse responses hh. . Ajn to estimate the data vector d using the extended 
algorithm. The received vectors d- • -Cm are processed by a channel equalizer 315 using 
the channel impulse responses h\..,hm to determine a spread vector S. A despreader 
320 using the transmission codes C despreads the spread vector S to estimate the data 
vector d. 

[0062] While the present invention has been described in terms of the preferred 
embodiment, other variations which are within the scope of the invention as outlined 

in the claims below will be apparent to those skilled in the art. 

* * * 
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